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Abstract 

The class of skewon-axion media can be defined in a simple and 
natural manner applying four-dimensional differential-form represen- 
tation of electromagnetic fields and media. It has been recently shown 
that an interface of a uniaxial skewon-axion medium acts as a DB 
boundary requiring vanishing normal components of the D and B vec- 
tors. In the present paper a more general skewon-axion medium is 
considered. It is shown that a planar interface of such a medium acts 
as a boundary generalizing both soft-and-hard (SH) and DB bound- 
ary conditions to SHDB conditions. Reflection of a plane wave from 
a planar SHDB boundary is studied. It is shown that for the two 
eigenpolarizations the boundary can be replaced by equivalent PEC 
or PMC boundaries. The theory is tested with a numerical example. 



1 Introduction 

The most general electromagnetic medium can be denned in terms of 36 
parameters, either in terms of four medium dyadics as 



D 9 
B 9 




(1) 
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in three-dimensional representation for Gibbsian vector fields denoted by D g 
, B 3 ,E 9 and H g , or in terms of a single medium dyadic as 

* = M|*, (2) 

applying four-dimensional differential-form representation for electromag- 
netic fields and media [1, 2]. In the latter case the field two-forms * and $ 
can be expressed as 

* = D - H A dr, $ = B + E A dr, (3) 

in terms of three-dimensional (spatial) two-forms D, B and one-forms H, E 
when d = ^£id Xi , i = 1..4 is the differential operator and r = ct is the 

normalized time. The medium dyadic M mapping two-forms to two-forms 
corresponds to a 6 x 6 matrix in any basis expansion of two-forms. The 
notation in this paper follows that given in [2]. 

The medium dyadic M has a natural (independent of any basis system) 
decomposition in three parts as 

M = Mi + M 2 + M 3 (4) 

and the dyadics Mj are respectively called principal, skewon and axion com- 
ponents of M [3]. The axion part is a multiple of the unit dyadic which in 
the notation of [2] can be expressed as 

M 3 = M 3 jW, (5) 

while both Mi and M 2 are trace-free dyadics. The components Mi and M 2 
can be defined so that the dyadics contracted by a quadrivector ejy = e 12 34 as 
e^Mi and e7v|_M 2 are respectively symmetric and antisymmetric. Physical 
properties of these components are discussed in [3]. 

The total number of 36 parameters is distributed by the three compo- 
nents so that the principal part Mi, corresponding to a trace-free symmetric 
6x6 matrix has 20 parameters, the skewon part M 2 , corresponding to an 
antisymmetric 6x6 matrix, has 15 parameters, and the axion part M 3 , has 
1 parameter. A medium consisting only of its axion parameter, M = M3 in 
(5), has been called PEMC, perfect electromagnetic conductor, because it is 
a generalization of both PMC (M 3 = 0) and PEC (1/M 3 = 0), [4, 5, 6]. The 
medium equations for the PEMC can be expressed in the simple form 

D = MB, H = —ME. (6) 
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2 Skewon-Axion Medium 

One can show that the most general medium with no principal component, 
the skewon-axion medium [8, 9, 10] (also called as the IB-medium [7]), can be 

expressed in terms of a dyadic B, corresponding to a 4 x 4 matrix involving 
only 16 parameters, as 

M = ( W = (W- (7) 

It will turn out that a plane wave in the most general skewon-axion medium 
is not governed by a dispersion equation, normally restricting the choice of 
the k vector in a given medium. Thus, basically, one can choose any k vector 
for a plane wave in such a medium. 

In [11] the problem of plane wave reflection from the planar interface of 
a uniaxial skewon-axion medium, defined by six medium parameters, was 
analyzed. It was shown that, despite the absence of any dispersion equation, 
the k vector in the skewon-axion medium was uniquely determined by the 
interface conditions. Moreover, it was shown that the interface conditions 
serve as boundary conditions 

n ■ D 9 = 0, n • B g = 0, (8) 

for Gibbsian fields at the interface of the uniaxial skewon-axion medium, 
dubbed as the DB conditions [12, 13]. It subsequently turned out that such 
conditions were defined already in 1959 [14] and that they are known to have 
unique mathematical properties [15, 16]. The DB conditions have been shown 
to play a central role in electromagnetic cloaking problems [17, 18, 19, 20]. 
More recently, the same DB conditions were shown to emerge at the interface 
of a simple skewon medium defined by just one parameter in the medium 
conditions [22] 

D = NB, H = NE, (9) 

which differ slightly those of the PEMC (axion medium) in (6). 

Fields in the skewon-axion medium can be compacly handled applying 
the four-dimensional formalism [2]. The Maxwell equations outside sources 
are 

dA$ = 0, dA* = 0, (10) 

whence the field two-form <fr can be (locally) represented in terms of a po- 
tential one-form <fi as 

$ = dA0. (11) 
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Substituting (2) and (7), the second Maxwell equation can now be expanded 
as 

dA* = dA(I^T T )|(dA</>) (12) 
= dA(I T |d) A</>) + dAdA(I T |(/>) (13) 
= dA((I T |d) A(f>) (14) 
= 0. (15) 

This shows us that the general solution can be represented as 

<f> = dip + W\d(p, (16) 

in terms of two scalar functions ip and ip, whence the field two-forms have 
the form 

$ = dA0 = dA(tl T |d)v?, (17) 
* = (I T |d) A</> = dA(I 2T |d)v?. (18) 

Instead of studying the general case, let us assume that the skewo-axion 
medium is defined by the dyadic 

1= B = \ + ae 3 + b(e 1 + Ae 4 ), (19) 

where £i,e 2 , s 3 are spatial one-forms and £ 4 = dr a temporal one-form, and 
together they form a basis. The corresponding reciprocal basis vectors are 
denoted by e, with the property e^e,- = 5ij [2]. 

Substituting (19) in (17) and (18), the field two-forms can be shown to 
satisfy the conditions 

£ 3 A (ei + As A ) A $ = 0, (20) 
e 3 A (ei + Ae A ) A * = 0, (21) 

for the special skewon-axion medium considered here. Inserting the expan- 
sions (3), the conditions become 

Ae 3 A e 4 A B + s 3 A £i A E A £ 4 = 0, (22) 
Ae 3 A £4 A D - £ 3 A £1 A H A £4 = 0, (23) 
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and they are equivalent with the spatial conditions 

Ae 3 A B + e 3 A £1 A E = 0, (24) 
Ae 3 A D — £ 3 A £1 A H = 0. (25) 

To conclude, any fields in a skewon-axion medium defined by (19) and (7) 
must satisfy the conditions (24) and (25). 



3 SHDB Boundary Conditions 

Let us interprete the result (24) and (25) in terms of three-dimensional Gibb- 
sian vectors applying the rules given in the Appendix. Denoting again the 
Gibbsian vectors by the subscript () g we can write 

e 123 |(£ 3 AB) = £ 3 |(e 123 LB) = e 3 • B„ (26) 

e 123 |(£ 3 AD) = £ 3 |(e 123 LD) = e 3 • D„ (27) 

e 123 1 (e 3 A £1 A E) = e 3 • (e x x E g ) = e 2 -E g , (28) 

e i23 1 (^3 A £1 A H) = e 3 • (ei x H fl ) = e 2 • H g . (29) 

Here we have applied the property e^Sj = Sij transformed to e, ■ = 5^ for 
the spatial basis vectors i,j = 1, 2, 3. 

Since the coefficient A in the conditions (24) and (25) has the dimension 
of velocity (due to the normalized time parameter r = ct), let us further 
express it as 

A = Tc=-^= = ^-, (30) 

in terms of a dimensionless scalar T. Here we have tacitly assumed that the 
fields are time harmonic as exp(jut) with k = u^fi e . Thus, the conditions 
(24) and (25) can be represented in the following form for the Gibbsian field 

vectors, 

Tooe 3 ■ B g + k Q e 2 ■ E g = 0, (31) 
Tcoe 3 ■ Ti g - k e 2 ■ H g = 0, (32) 

and they are satisfied by any fields in the skewon-axion medium under consid- 
eration. One may note that the pair of conditions is invariant in the duality 
transformation D g B 9 , E 9 -h- — H 5 which leave the Maxwell equations 
invariant. 



6 



Assuming now a planar interface at e 3 • r = 0, with normal unit vector 
e3, between the skewon-axion medium and an isotropic medium (parame- 
ters n ,e ), the conditions (24) and (25), and their Gibbsian counterparts 
(31) and (32), are actually continuous across the interface. Since these con- 
ditions at the interface are enough to determine the fields in the isotropic 
medium without having to solve for the fields behind the interface, they can 
be considered as boundary conditions. Although (31) and (32) were intro- 
duced through a consideration of a planar interface, their form is applicable 
to curved boundaries as well, with e 3 and e 2 denoting respective unit vectors 
normal and tangential to the boundary surface. 

The boundary defined by (31), (32) appears as a generalization of both 
the DB boundary defined by the conditions (8) and the soft-and-hard (SH) 
boundary [23, 24, 20] defined by the conditions 

e 2 • E g = 0, e 2 • H 9 = 0. (33) 

In fact, the SH conditions are obtained for the limiting parameter value 
T — > while the DB conditions are obtained from (31), (32) for T — > oo. 
For convenience, the more general boundary defined by (31) and (32) and 
combining both SH and DB conditions will be called the SHDB boundary. 

One can immediately observe that, due to its self-dual property, an object 
with rotational symmetry and SHDB boundary should appear invisible to the 
radar because the back-scattering cross section of such an object is zero [21]. 

4 Gibbsian Medium Conditions 

The skewon-axion medium was defined in terms of four-dimensional quan- 
tities through (7), (19). For the convenience of readers not familiar with 
the four-dimensional formalism and for double-checking the results, let us 
now find the relation between the skewon medium interface and the bound- 
ary consitions (31), (32) in terms of three-dimensional Gibbsian vectors and 
dyadics. For that purpose we first find the Gibbsian medium conditions by 
expanding 

(I^T) T = 2B = \ {2)T - (e 3 A F A a + (ei + Ae 4 ) A F A b), (34) 
and expand 



(e 3 AF Aa)|$ = (e 3 AF Aa)|(B + EAe 4 ) 
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= e 3 A (ajB) + (e 3 A E)(a|e 4 ) - (e 3 A e 4 )(a|E). (35) 

The expansion for the second term is obtained similarly, replacing e 3 by 
£\ + Ae 4 and a by b. Let us also expand the two vectors in their spatial and 
temporal components as 

a = a s + a 4 e 4 , b = b s + 6464, (36) 

With these we can split the medium equation 

D-HAe 4 = (IaT) t |(B + EA£ 4 ) (37) 
in its spatial and temporal parts and identify the fields as 

D = 2BB — £ 3 A (a s J B) - e 1 A {b s \ B) 

-(a 4 £ 3 + & 4 £i)AE, (38) 

H = -Ab s \B-{Ab 4 + 2B)E 

-(s 3 a s + £ib s )|E. (39) 

Applying the transformation rules given in the Appendix, the medium 
conditions can be written for the Gibbsian vectors as 

B g = 2BB g - e 3 x (a s x B fl ) - e l x (b s x B g ) 

-(a 4 e 3 + 6 4 ei) x E g , (40) 

H g = -Ab s xB g - (Ab 4 + 2B)~E g 

-(e 3 a s + eib a ) • E r (41) 

Expressing the Gibbsian medium conditions in the form 



where the medium dyadics are understood in the Gibbsian sense, their ex- 
pressions become 

!! = 2B\ g - e 3 x (a s x \ g ) - ei x (b s x \ g ), (43) 

e = -(a 4 e 3 + b 4 e 1 ) x \ g , (44) 

W 1 = -Ah s x = \ g , (45) 

I = -(Ab 4 + 2B) = \g- (e 3 a s + ei b s ). (46) 
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One must notice that the dyadics e and JT 1 are antisymmetric and, conse- 
quently, do not have inverses. Thus, it is not possible to express the skewon- 
axion medium equations in the form (1). 

The most general skewon-axion medium can be defined in terms of three- 
dimensional Gibbsian vectors and dyadics in the form [7, 11] 

D 9 = trA B 9 -A-B 5 -cx E g , (47) 
H 9 = -gxB 9 -F-E 9 -aE 9 , (48) 

where A is a dyadic, c, g are two vectors and a is a scalar. Together they 
make 9 + 3+3 + 1 = 16 scalar parameters. Comparing with the representation 
(43) - (46), we can identify the quantities in (47) and (48) as 

t = B\ + a s e 3 + b sei , (49) 

c = 6 4 ei + a 4 e 3 , (50) 

g = ^4b s , (51) 

a = B + Ah, (52) 

which shows us that the number of free parameters of the medium defined 
by (19) must be less than 16. Expressing the dyadic A as 

t = Be 2 e 2 + (Be 3 + a s )e 3 + (Be 1 + b s )ei, (53) 

its definition is seen to require 1 + 3 + 3 = 7 parameters in a given basis 
ej, while c and g require 2 and 1 parameters, respectively. Thus, the total 
number of free parameters in the medium under consideration appears to be 
10. The pure axion (PEMC) medium (6) special case is obtained for 

= M- 

A=yl„ c = 0, g = 0, a = M/2. (54) 

5 Reflection from SHDB Boundary 

Let us consider the basic problem, reflection of a plane wave from the planar 
interface of the skewon-axion half space z < 0, in terms of Gibbsian quantities 
corresponding to time-harmonic fields. 
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5.1 Dispersion Equation 

Denoting for simplicity 

p = k/w, (55) 
a plane wave with the dependence exp(— juup ■ r) satisfies 

pxE s = B s , p x H 9 = -D r (56) 

To obtain a condition for a plane wave in the skewon-axion medium, substi- 
tuting (47) and (48) and eliminating B g in (56) leaves us with an equation 
of the form [7] 

U(p) • E g = 0, (57) 

where we denote 

^(p) = q(p) x I, (58) 
q(p) = (g • P - a)p - c + p • t. (59) 

Because detD(p) = for any p, there is no dispersion equation to restrict the 
choice of p. Thus, a plane wave with any vector k is possible in the skewon- 
axion medium. One can further show that q(p) = for any p exactly when 
the medium has no skewon component, i.e., when it is a pure axion (PEMC) 

medium with c = g = and A = a\ s . In the following we exlude the pure 
axion medium from the analysis. 

5.2 Field Conditions 

The field vectors E g and B g of the plane wave can be expressed as 

E 9 = q(p) = (gp-a)p-c + pl, (60) 
B g = p x E g = -p x c - p • A x p. (61) 

Inserting (49) - (52), the field expressions for the special skewon-axion medium 
become 



Eg = (b s • p - b A ) (Ap + ei) + (a s • p - a 4 )e 3 , (62) 
B g = (b s • p - 6 4 )p x ei + (a s • p - a 4 )p x e 3 . (63) 
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The form of the components 

E g -e 2 = A(b a -p-6 4 )(p-e 2 ), (64) 
B 9 -e 3 = -(b s • p - 6 4 ) (p • e 2 ), (65) 

implies that, for the general vector p, the fields satisfy the condition 

Ae 3 ■ B g + e 2 • E g = 0, (66) 

which coincides with (31) when (30) is taken into account. If p satisfies 
p ■ e 2 = 0, we obtain both B s ■ e 3 = and E 9 ■ e 2 = 0, which are actually the 
same condition because they are connected by 

e 3 • (p x E g - B g ) = (ei • p)(e 2 • E g ) - e 3 • B g = 0. (67) 

Inserting (60) and (61) in (47) and (48), with (49) - (52) and after some 
algebraic steps, the following expansions for the other field components are 
obtained, 

e 3 -D s = -(p • e 2 )((a s • p - a 4 )(b s ■ e 3 ) + 

+(b s p- 64) {Ab 4 + 2B + b s -e 1 )) (68) 

e 2 -H s = -A(p ■ e 2 )((a s • p - a 4 )(b s • e 3 ) + 

+(b s -p-b 4 )(Ab 4 + 2B + b s -e 1 )). (69) 

Thus, in the general case, the fields satisfy 

Ae 3 ■ E> g - e 2 • H g = 0, (70) 

which coincides with (32) when (30) is taken into account. Again, for p e 2 = 
we obtain e 3 • D 3 = and e 2 • H 9 = 0, which are the same condition. 

It now follows that, inserting (30), the sum of incident and reflected plane- 
wave fields at the boundary of the isotropic medium satisfy the conditions 



r W (Bj + BJ) • es = -k (Ei + E r g )-e 2 , (71) 

7MD*+Dp.e 3 = fc (H*+Hp.e 2 , (72) 

which can be rewritten as 

T(k l xE; + k r xE;).e 3 = -A; (E; + E;)-e 2 , (73) 

T(k l x H; + k r x Hp • e 3 = -A; (H; + Hp-e 2 . (74) 



11 



Applying the relations 

Hj = —Y x Ej, H; = — k^ x E r g , (75) 

k r = (eiei + e 2 e 2 - e 3 e 3 ) • k\ (76) 

k l • E; = 0, k r • E^ = 0, (77) 

k r • k r = k l • k* = k 2 Q , (78) 

one can find relations between the fields E^ and E^ which depend on the 
medium parameters and the vector k* of the incident wave. Skipping the 
algebraic details, we obtain the following result for the field components 
parallel to the boundary: 

(79) 



2 / V ^2 



with matrices defined by 



a = Tk\ (81) 

p = k + Tk[ (82) 

7 = k 2 +Tk k\-kf (83) 

5 = (k{ + Tk )ki (84) 

The determinants of the two matrices become 

A = detA = -detB 

= ko ((Tk + k\) 2 -(T 2 -l)kf). (85) 

The reflected field components can be solved as 

E[ = -El + ^Ei + SEi) (86) 

El = -E^—^El + SEl) (87) 

El = Ei + ^aki + Pki^El + SEi) (88) 
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As a check, for T = we obtain 

2Af 

E l = ~ E 2, E i = E\+ i2 J_ l i2 E l 2 . (89) 

These coincide with the SH conditions because the latter relation equals 
e 2 • = — e 2 • Hg, as can be verified. Also, for T — >■ oo we obtain Eg — > —£3 
and ifg — >■ — if 3 which correspond to the respective DB conditions e 3 -D 9 — > 
and e 3 • B 9 — > for the total fields. 

As another check we may consider the normally incident plane wave with 
k* = — k Q e 3 . In this case we obtain E{ = E\ and E r 2 = —E l 2l which coincide 
with the SH conditions. This is independent of any finite value of the pa- 
rameter T. Thus, if 1/T is a small quantity, the SHDB boundary acts as a 
DB boundary, except for normal incidence where the SH condition suddenly 
takes over. 

Finally, for a plane wave satisfying e 2 • k* = both SH and DB conditions 
are simultaneously satified for any T and k\ values. 



5.3 Eigenfields 

To find the eigenpolarizations satisfying 

E[ = XE{, E r 2 = XEi (90) 

we must solve 

det(A - \B) = 0, (91) 
for A. The solutions are simply 

A± = ±1, (92) 

whence the eigenfields satisfy either PMC (A + ) or PEC (A_) conditions at 
the SHDB boundary. The eigenfields satisfy 

ai E\ + = (3Ei + , (93) 
7 £i_ = -8Ei_. (94) 

One can show that these equal the conditions 

Twe 3 • B; + + k e 2 ■ E; + = (95) 
Tcue 3 ■ D;_ - k e 2 ■ IT. = 0, (96) 



13 



which restrict the polarization of the two eigenwaves. Considering the two 
special cases, for T = (SH boundary) the eigenfields are either TE or TM 
with respect to e 2 while for T — > oo (DB boundary) they are TE or TM with 
respect to e 3 . Actually, in an isotropic medium, the polarization conditions 
are of the form 



which can be understood as generalizations of TE/TM conditions. Such 
polarization conditions have been previously encountered in media known by 
the name decomposable media, in which any fields can be decomposed in two 
non-interacting components [25]. 

Since (95) and (96) are also valid for the respective reflected eigenwaves, 
they are valid for the total eigenfields. The conditions are linear in the field 
vectors and independent of the k vector of the plane wave. Thus, they are 
valid for linear combinations of plane waves and, actually, for any fields. 
Thus, any field can be decomposed in two eigenfields satisfying either (95) 
or (96) for which the SHDB boundary can be replaced by respective PMC 
and PEC boundaries. 

To show that any given field can be expanded as a sum of eigenfields, let 
us consider the general plane wave and express the conditions (95) and (96) 
as 



a+ ■ E g+ + b+ ■ H 9+ = 
a • E 9 _ + b • H y = 0, 



(97) 
(98) 



Te 3 • k x E 9+ + k a e 2 • E 9+ 
Tk Q e 3 • E 9 _ - e 2 • (k x E 9 _) 



= 



(99) 
(100) 



which are of the form 



c 3 -E, 



0, c 2 • E,_ = 0, 



(101) 



with 



c 3 = Te 3 x k + k D e 2 , 
c 2 = Tk Q e 3 - e 2 x k. 



(102) 
(103) 



These vectors satisfy 

e 3 • c 3 = 0, e 2 • c 2 = 0. 



(104) 
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Expanding a = k x ((c 3 x c 2 ) x E s ) in two ways as 

a = k x c 2 (c 3 • E g ) - k x c 3 (c 2 • E g )) (105) 
= (c 3 xc 2 )(k-E fl ) -(k-(c 3 xc 2 ))E s , (106) 

and applying k • E g = 0, the eigenfield representation 

E 9 = E fl+ + E 5 _, (107) 

is obtained by denning 

„ k x c 3 c 2 „ 

= M^t)' E - < 108 > 

= k»cc, 
k • (c 3 x c 2 ) 

The expansion requires k ■ (c 3 x c 2 ) ^ which corresponds to 

T^kl + kD + kf + kl^O. (110) 

For the magnetic field H g = H 3+ + H 9 _, the polarization conditions can 
be shown to take the form 

c 2 -H 3+ = 0, c 3 -H s _ = 0, (111) 

whence the magnetic eigenfields obey the expressions 

Obviously, the two eigenfields are dual to each other because the transfor- 
mation E g — > H g , Hg — > Eg corresponds to E g+ — > E 9 _, E 9 _ — > E g+ . 



6 Numerical example 

To have an idea of SHDB boundary in action, let us consider as an example 
a plane wave incident to the SHDB plane with (p = n/2 in Figure 1. Thus, 
the incident wave vector satisfies 



ei • k l = 0, k l = &; o (e 2 sin0 + e 3 cos#). 



(114) 
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Figure 1: Geometry of the reflection problem. A plane wave is incident on the 
SHDB boundary with directions defined by the angles 6, (p of the spherical 
coordinate system. The fields satisfy k* • = and k r • = 0. 



The relation between the incident and reflected field components transverse 
to e 3 can be represented by the planar reflection dyadic R as 

e 1 E[ + e 2 E r 2 =^- ( ei E{ + e 2 Ei), (115) 

R = eieii?n + eie 2 -Ri2 + e 2 eifl 2i + e 2 e 2 fl 22 , (116) 



with 



cos 2 9 - T 2 sin 2 9 

cos 2 9 + T 2 sin 2 6'' 
2Tsm6 

cos 2 6 + T 2 sin 2 #' 

2Tsin#cos 2 fl 
cos 2 # + T 2 sin 2 #' 
fl 22 = -R11. (120) 



„ J_ Dill V 

«« = -ynr-^r-TT,. ("8) 



The reflection dyadic satisfies 

trl = flu + fl 22 = 0, (121) 
detl = -R 2 U - fli 2 fl 2 i = -1, (122) 
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as is also obvious from the eigenvalues A± = ±1. 

Dependence of each component Rij of the reflection dyadic on the angle of 
incidence 6 is depicted in Figure 2 for three parameter values T = (0.1, 1, 10). 
According to the previous theory, for T = 0.1 the SHDB boundary should act 
almost like the DB boundary while for T = 10 it should rather resemble the 
SH boundary. This is clearly seen from the graphs of Rn and R 22 because 
for T = 0.1 (solid curves) the component E 2 reflects almost from PEC with 
R 22 ~ — 1 and the component E\ almost from PMC with Rn ps +1. There 
is a deviation for angles close to grazing. The opposite is the case for T = 10 
with deviation close to normal incidence which is a known effect for the DB 
boundary. There are small cross-polarized reflections for small and large 
values of T. 

For T = 1 the SHDB boundary exhbits properties not shared by SH or 
DB boundaries. In this case we have 

R u = -R 22 = cos 26, R 12 = 2 sin 6, R 21 = sin 26 cos 6. (123) 

For example, for the incidence angle 6 = ir/4 the reflected field appears 
totally cross polarized, Ru = R 22 = 0. This is due to the nonreciprocal 
character of the SHDB boundary, which stems from the gyrotropic medium 
dyadics of the skewon-axion medium, (43) - (46). 

As an another example we could consider the case e 2 ■ k* = 0, or <p = in 
Figure 1. Forming the reflection dyadic, we would obtain R n = —R 22 = 1 
and R\ 2 = R 2 \ = for all 6 and T values. This would make a dull diagram. 
As was already pointed out, the case e 2 ■ k l = corresponds to reflection 
from either an SH boundary or a DB boundary, both of which yield the same 
result. Actually, the DB boundary can be characterized as an isotropic soft 
surface [11], while the SH boundary is a soft surface for waves arriving from 
a certain direction, which in the present case corresponds to e 2 ■ k l = 0. 



7 Conclusion and Discussion 

The class of skewon-axion media can be described mathematically in natu- 
ral (coordinate-independent) manner applying four- dimensional formalism. 
Instead of the 36 medium parameters of the most general linear medium, 
represented by a 6 x 6 matrix, the skewon-axion medium only involves 16 
parameters defined by its characteristic 4x4 matrix. General solutions to the 
Maxwell equations can be expressed in simple form in the four-dimensional 
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Figure 2: Components of the reflection dyadic Rij for a plane wave incident as 
ip = 7r/2. The angle of incidence 9 is shown in degrees. Curves in each figure 
correspond to parameter values T = 0.1 (solid blue line); T = 1 (long-dashed 
red line); and T = 10 (short-dashed green line). Note that the co-polarized 
reflection coefficients satisfy R u = — in- 
form. Considering a plane wave in terms of three-dimensional Gibbsian vec- 
tor fields, it is shown that there is no characteristic equation for the vector 
k in the medium which leaves a lot of freedom for the definition of the wave. 

In this paper a certain class of skewon-axion media is studied whose num- 
ber of parameters is limited to 10. Such a medium has three specified axes 
denoted by the constant orthogonal unit vectors ei,e2 and e 3 . It is shown 
that, in such a medium, two scalar conditions (31), (32) are satisfied by 
the four Gibbsian field vectors. The form of the conditions suggests that, if 
there is a planar interface orthogonal to the vector e3, (31) and (32) serve as 
boundary conditions for the fields. Since the conditions generalize those ob- 
tained for the soft-and-hard (SH) boundary and for the DB boundary for the 
limiting cases of the free parameter T, the novel boundary has been dubbed 
as the SHDB boundary. 

Plane wave reflection from a SHDB boundary was considered and ana- 
lytical expressions for the reflected field components were derived. It was 
shown that, for two eigenpolarizations, the SHDB boundary can be replaced 
by effective PEC and PMC boundaries just like for the SH and DB bound- 
aries. The theory was tested by a numerical example demonstrating the 
nonreciprocal property of the SHDB boundary. 
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Although it appears quite easy to define the generalized boundary condi- 
tions (31) and (32), quite simple in form, one may note that their realization 
in terms of a medium interface is not so obvious. In fact, it appears that, 
without the four-dimensional formalism, the introduction of the skewon-axion 
medium is a difficult task. Also, it is not at all obvious that a medium de- 
fined by (47) and (48) with (49) -(52) will yield a set of boundary conditions 
since it does not pop out directly from the medium equations but requires 
consideration of fields in the medium. 

Realization of the skewon-axion medium by some metamaterial is a chal- 
lenge for which there may not exist an easy solution. However, since the 
boundary conditions (31) and (32) appear more important than the medium 
itself, it should be easier to realize the SHDB boundary, whenever a useful 
engineering application for the boundary has been discovered. A correspond- 
ing realization for the DB boundary was recently found in terms of a planar 
structure containing metamaterial inclusions [26]. 

Appendix: Gibbsian Representation of Forms 

The connection between spatial multiforms and Gibbsian vectors can be ex- 
pressed according to the following rules which are based on the spatial metric 
dyadic 

\ g = eiei + e 2 e 2 + e 3 e 3 , (124) 

which serves as the Gibbsian unit dyadic. Thus, a spatial one-form a s = 
a.\£\ + a 2 e 2 + «3^3 is transformed to a Gibbsian vector as 

ct s ->■ ct g = \g\a s = aid + a 2 e 2 + a 3 e 3 . (125) 

Similarly, a spatial two-form B = B 12 e 12 + B 23 e 23 + B 31 e 31 is transformed as 

B ->■ B g = e 123 [B = B 12 e 3 + B 23ei + B 31 e 2 . (126) 

The bar product of a spatial vector a s and a spatial one-form ct s transforms 
to a dot-product of Gibbsian vectors, 

a s |ai s ->• a s ■ a g = a s -\g\ais = a. s \a s . (127) 

The wedge product of two spatial one- forms ct s , (3 S becomes the cross product 
of two Gibbsian vectors, 

a s A(3 s ->■ a g x f3 g = e 123 [(a s Af3 g ), (128) 
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and the wedge product of three one-forms transforms to a scalar, 

<x s A(3 s A ls -> <x g -(f3 g x lg ) = as \(e 123 l((3 s A ls )). (129) 

Finally, the contraction product of a spatial vector a s and a spatial two-form 
B is transformed to the cross product of the Gibbsian vectors a s and B 9 as 

ajB -> a s x B 9 =l 9 |(a s jB). (130) 
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